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ALMOST SURE CENTRAL LIMIT THEOREM FOR BRANCHING 
RANDOM WALKS IN RANDOM ENVIRONMENT 

By Makoto Nakashima 

Kyoto University 

We consider the branching random walks in d-dimensional in- 
teger lattice with time-space i.i.d. offspring distributions. Then the 
normalization of the total population is a nonnegative martingale and 
it almost surely converges to a certain random variable. When d > 3 
and the fluctuation of environment satisfies a certain uniform square 
integrability then it is nondegenerate and we prove a central limit 
theorem for the density of the population in terms of almost sure 
convergence. 

1. Introduction. We write N = {0, 1, 2, . . . , }, N* = {1, 2, . . . , } and Z = 

{±x:x E N}. For x = (xi, . . . , rr^) S M*^, |x| stands for the i^-norm: \x\ = 

Eti kd- For C = i^,),^^, G R^', 1^1 = Z.ez" \U Let {n,T,P) be a prob- 
ability space. We write E[X] = f X dP and E[X ■.A\ = f^X dP for a random 
variable X and an event A. We denote the constants by C,Ci. 

We consider the branching random walks in random environment. Branch- 
ing random walks have been much studied [1, 2] and a central limit theorem 
for the density of the population has been proved in the nonrandom environ- 
ment case [2]. Also, in the random environment case, one has been proved 
in the sense of "convergence in probability" [20] when d>3 and the fluctu- 
ation of environment is well moderated by the random walk. In this article 
we prove a central limit theorem in the sense of "almost sure convergence" 
under the same condition as in [20] . The time-space continuous counterpart 
is the branching Brownian motion in random environment for which the cen- 
tral limit theorem has been proved in [15]. On the other hand, a localization 
property has been proved in [10] for the branching random walks in random 
environment if the randomness of the environment dominates. 
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It has been mentioned that the branching random walks in random en- 
vironment (BRWRE) have a similar structure to the directed polymers in 
random environment (DPRE) [3, 5, 7, 20]. Also, we will see the relation 
between BRWRE and DPRE in Section 1.3. A central limit theorem has 
been proved for a Markov-chain-generalization of the directed polymers in 
random environment [4, 6, 12, 13] assuming a certain square integrability. 
Since we use an analogy to [13], we extend the framework to contain the 
branching random walks in random environment. 

1.1. Branching random walks in random environment. We consider par- 
ticles in Z'^, performing random walks and branching into independent copies 
at each step of the random walk: 

(i) At time t = there is one particle at the origin x = 0. 

(ii) When a particle is located at site x G Z'^ at time t E N, it moves to a 
uniformly chosen nearest neighbor site and is replaced at time i + 1 by 
fe-particles with probability qt,x{k){k E N), 

where we assume that the offspring distributions qt^x = {qt,x{k))k&n are i.i.d. 
in time-space {t,x). This model is investigated in [3] and we call it the 
branching random walks in random environment (BRWRE). Let Nt^y be 
the number of the particles which occupy the site y E Z*^ at time t. Let Nt 
be the total population at time t. In this article we study the behavior of the 
density pt{y) = l{/vrj>o}- We look at the branching process to give a more 
precise definition of the branching random walks in random environment. 
First, we define 14, n E N, Vm by 

Vo = {l}, Vi = {N*f,...,Vn = {n*r~^' forn>l, 
Vn = U Vn. 

nGN 

Then we label all particles as follows: 

(i) At time t = there exists just one particle which we call 1 E Vq. 

(ii) A particle which lives at time t is identified with a genealogical chart 
y = (1, yi , . . . , yt) E . If the particle y gives birth to ky particles at time 
t, then the children are labeled by (1, yi, . . . , yt, 1), . . . , (1, yi, ■ ■ ■ ,yt, ky) E 
Vt+i. 

By using this naming procedure we rigorously define the branching random 
walks in random environment. This definition is based on the one in [20]. 
Note that the particle with name x can be located at x anywhere in Z'^. 
As both information genealogy and place are usually necessary together, it 
is convenient to combine them to is^= (x,x); think of x and x written very 
closely together. 
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• Spatial motion. A particle at time-space location (t,x) is supposed to 
jump to some other location {t + l,y) and is replaced there by its children. 
Therefore, the spatial motion should be described by assigning a destina- 
tion for each particle at each time-space location {t,x). So we are guided 
to the following definition. Let the measurable space {^Ix^J'x) be the set 
^^d^NxZ'^xVn with the product a-field and Qx B X = (^t^^)(t,^,y)gNxZ<ixVN- 
We define Px G V{Qx ,J^x) as the product measure such that 

PxiXl = e) = lM' ^^l"l = l' 
lo, if|e|/l 

for e G Z'^ and {t,x,y) G N x Z"^ x Vn- Here we interpret Xf^ as the step at 
time t + 1 of the particle y at time-space location {t,x). 

• Offspring distribution. We set Qq = 7^(N)^^^'* where V(N) denotes the 
set of probability measure of N: 

P(N) = L= iqik))km G [0, lf;Yl l^k) = 1 

km 

Thus each G is a function {t,x) i— )• qt^x = iQt,x{k))km from N x Z'^ to 
V(N). We interpret qt^x as the offspring distribution for every particle oc- 
cupying the time-space location {t,x). The set 'P(N) is equipped with the 
natural Borel cj-field introduced from that of [0,1]^. We denote by J-q the 
product (T-field on ilq. 

We define the measurable space {nK,J^K) as the set N^^^ with the 
product (7-field and Q,k B K = (-f^t^2:)(f,a;,y)GNxZ'*xVN- each fixed q £ ^Iq 
we define P^ G V{^k,^k) as the product measure such that 

p9 (i^y^ = k) = qt,x{k) for all {t, x, y) G N x Z'^ x Vn and A: G N. 

We interpret Kf^ as the number of the children born from the particle y at 
time-space location {t,x). 

We now define the branching random walks in random environment. We 
fix a product measure Q G V{VLq,Fq), which describes the i.i.d. offspring 
distribution assigned to each time-space location. In the following we also 
use Q as Q-expectation, that is, we write Q\Y] = jYdQ and (5[y:y4] = 
jj^YdQ for a g-random variable Y and an Jg-measurable set A. Finally, we 
define by 

= Vlx X VLk X 

and P'i,PeV{n, F) for g G by 
P'' = Px®P'k®5q, 



F = Fx® Fk ® Fq 



I Q{dq)P''. 
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We want to look at Nt^y but here we investigate more detailed information. 
We define Nj^ by 

(1.1) N^^ = Ijthe particle y is located at time-space location 

for (t,x,y) G N X Z'^ X Vn- Here we set N^^ = 5q'^ where 6 is the Dirac 
function such that 

^i.y^fl, if x = 0,y = 1 G Vb, 

'''^ 1 0, otherwise. 

Then we can describe Nj^ inductively by 

(1.2) = N^-i,.Hy - ^ = xti,., 1 < y/x < Kti,.} for t > 1, 

xGVn 

where y/x is given for x,y G Vn as follows: 

f k, if X= {l,Xi,. . . ,Xn) £Vn,y = {l,Xi,. . .,Xn,k) £Vn+l 

y/x = < for some n G N, 

I oo, otherwise. 

Moreover, Nt^y and Nt can be rewritten respectively as 

(1.3) ^t,y=YNly and iV* = J] A^^, 

yeVN 

for t G N,y G Z"^. We remark that the total population is exactly the clas- 
sical Galton-Watson process if qt^x = Q, where q G 'P(N) is nonrandom. For 
simplicity we write (1.2) as 

(1.4) = ^ N^.,,,Al'ly for t > 1, 

xGVpj 

where we set 

Atly = i{y-x = xti,x, 1 < y/x < J. 

This formula is similar to the one for directed polymers in random environ- 
ment and linear stochastic evolutions. For p > 0, we write 

m(P) = Q [mg] with = ^ Fg*,^ (A;) , 

fceN 

m = m^^^ , mt^x = '"T'f x • 

We set 

(1.5) Nly = Nly/m\ Nt,y=Nt,y/m' and Nt = Nt/m' 
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for [t, y, y) G N x Z'^ x Vn- We prove later that A^^ is a martingale with respect 
to Tt = cr{As ■.s<t) where At = {A^'^y : {x, x), (y, y) G N x Vm}- Therefore, 
the following limit always exists (see Theorem 1.2): 



A^oo = lim Nt, P-a.s. 

t~+oo 



It is easy to see that 
(1.6) 



h^lU)^ if = l,y/x = A;,A;G: 



2d^ 

0, otherwise, 



and 



a^^^ = m for X G Z'^, X G Vn, 

where q{j) is the Q-expectation of qt,x{j)- 

1.2. Properties. In this section we look through the properties of BR- 
WRE. First we introduce an important Markov chain and represent Njy by 
using it. We define the Markov chain (5,Pf) = {{S^ ^S'^),?^^^^) on x Vn 
for X = (x,x) G Z'^ X Vn, independent of {At}t>i^ by 

P|(5o = ^) = P|;>(So = (x,x)) = 1, (x,x) G Z'^ X Vn 

and for each x, y G U^, x, y G Vn, 

Ps{St+i = {y,Y)\St = {x,^)) 

M m 

0, otherwise. 

We remark that we can regard and S'^ as independent Markov chains 
on Z'^ and Vn, respectively, and that is a simple random walk on U^. 
Here we introduce a certain martingale which is essential to the proof of our 
results: 

jiSs 

Co = l and fort>l Ct = H 



l<s<t 



s-1 



9 5^ S^/S^ 

where A^*g^__^ = A^*^j^' and 05^/5^ .^ ^^^i.^r^ • In fact, this is a martin- 
gale with respect to the filtration defined by Tit = cr{Au, Su, u<t) as in the 
following lemma where for t = 0, Hq = (t(5o). 
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Lemma 1.1. Q is a martingale with respect to %t- Moreover, we have 
that 

(1.7) Nly = m'El[(:t;St = {v,y)], P-a.s. for t en,{y,y) eZ^ x Vn, 
where Eg[-] denotes the expectation with respect to Pg'= P^'^. 

Proof. Since {At}t>i are i.i.d. random variables, it follows from the 
independence of {At}t>i and {S'f}t>i that 



aSi 



t~l, 



where ^ is the product probability measure of P and Pg and where 
E\ g denotes the expectation with respect to P^g- We now prove (1.7) by 
induction. It is easy to see that (1.7) holds for t = 0. If (1.7) holds for t > 0, 
then 



x,y 



Since we have that 



A 



t+l,x,y 



mEr 



aSi 



\Si = {y,y) 



«5i/5o 

(1.7) holds for t + 1 and the proof is complete. □ 



We remark that Nt^y = m}E^[C,t]Sl = y\. From this lemma we obtain an 
important result. The following theorem means that a phase transition oc- 
curs for the growth rate of the total population. 

Theorem 1.2. Nt is a martingale with respect to Tt = a{As : s <t) and 
there exists the limit 



N. 



lim Nt 



OO lliiJ. -IV 4, 



P-a.s., 



(1.8) 

and 

(1.9) E\N^] = l or 0. 

Moreover, E\Ncc\ = ^ if o,nd only if the limit (1-8) is convergent in L^{P). 
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Before we prove Theorem 1.2 we introduce some notations and definitions. 

For (s, z) e N X Z'^ X Vn, we define N,''''^ = (A^t^;,'y')fe,y)gz^xVM and W,'''^ = 
'y'y''{!/.y)GZ'*xVN' * ^ respectively by 

^s,2,z ^ ey,z ^ f 1, if y = z, and y = z, 
o,y,y |o, otherwise, 

(1.10) iVS,y= E and 

We remark that we can regard iV*'^'^ = {-/Vs,'y,y}(y,y)gz<ixVN as the state of 
the branching random walks starting from particle z at time-space {s,z) 
observed at time s + t. 

Proof of Theorem 1.2. The limit (1.8) exists by the martingale con- 

vergence theorem since Nt is a nonnegative martingale and i = E[Noo] < 1 
by Fatou's lemma. To show (1.9) we will prove that i = i"^ using the argu- 
ment in [11]. With the notation (1.10) we write 

where A^^ ^ is defined by (1.5) for (s, z, z) G N x Z*^ x Vn- Since every (s, z, z) G 
N X Z"^ X Vn, N^'^'^ is a martingale with respect to J-'f = a{As^u ■.u<i) and 
has the same distribution as Nt, the limit 



iv:;^'^= lim Nr" 

exists almost surely and is identically distributed as A^oo- Moreover, by let- 
ting t — )• oo, we have that 

-'^oo / ^ s,z-^^ oo 

^GZ^.zGVn 

and hence, by Jensen's inequality, that 

E[eic^{-N^)\Fs\ > e-x.v{-E\N^\Fs\) = exp(-iV,£) > exp(-¥,). 
By letting s — )■ oo in the above inequality, we obtain 

exp(-iVoo) >' exp(-iVoo^) > exp(-iVoo) 



and thus A^oo A^oo^- By taking expectation we get i = . Once we know 
(1.9), the final statement of the theorem is standard (e.g., [9], formula (5.2), 
pages 257-258). □ 
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We refer to the case E[Noo] = 1 as the regular growth phase and to the one 
-E'[-^oo] = as the slow growth phase. The regular growth phase means that 
the growth rate of the total population is the same order as its expectation 
m* and the slow growth phase means that, almost surely, the growth rate is 
slower than the growth rate of its expectation. 

We discuss the case of the regular growth phase in this article. The slow 
growth phase is partially studied in [10]. If Nt is uniformly square integrable 
then it is the regular growth phase since Nt is a martingale. 

Here we give the main theorem in this article. 

Theorem 1.3. Suppose that d>3 and 

(1.11) m>l, m(2)<oo and ^ 1^ 

where tt^ is the return probability of a simple random walk in Z*^. Then for 
all f £Cb(M.'^), 

(1.12) hm V /f4l^M = ^oo / f{x)duix), P-a.s., 

where Cb{M.'^) stands for the set of bounded continuous functions on R'^ and 
V is the Gaussian measure with mean and covariance matrix 

The proof of Theorem 1.3 will be given in the next section. 

Remark. From Lemma 1.1 we can rewrite (1.12) as 



(1.13) lim Pg 



t—>oo 



Noo / f{x)diy{x), P-a.s. 



In Lemma 2.2 we see that (1.11) is equivalent to sup^^i E[{Nt)'^] < oo so 
that we have E[Noo] = 1, that is, P{Noo > 0) > 0. Also, if we set pt{x) = 
^^l{Nt > 0}, we can interpret pt{x) as the density of the particles. From 
this observation we can regard Theorem 1.3 as a central limit theorem for 
probability measures with the density pt{\/tx) on {A^oo > 0}. 

1.3. Relation to directed polymers in random environment. In the end of 
this section we discuss the relation between BRWRE and DPRE (see [20], 
pages 1631-1634, for more detailed information). 

• Random walk. [S[,P'^i) is a simple random walk on d-dimensional lattice 
defined on the canonical path space (Os-/, -Fg'). Pg, is the unique probability 
measure on {VLs',J^s') such that S'^ — Sq, . . . , 5j — S[_^ are independent and 

P|,(5^ = x) = l, Pl,{Si_,,-S't = e)-- 
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where e G Z'^. For x = we write simply Ps' as Pg,. We denote by Eg, 
the P|,-expectation. We can regard {S'^^Pg,) as an independent copy of 

• Random environment, rj = {rit,x : G N x Z*^} are M- valued i.i.d. ran- 
dom variables which are nonconstant and defined on a probability space 
{VLr^,Fri,Q') such that 

gA{/3) drf Q'[exp(/37?t,^)] < oo for ah /? E R. 

• Polymer measure. For any t G N, define the polymer measure on 
{^S',^S') by 

diit = ^fiMl^Ht-t\{P))dPs', 

where /3 > is a parameter and 
t-i 

Ht = Y,^u,s'^ aiid Zt = E%[eMm-tKm 

u=0 

are the Hamiltonian and the partition functions. We call this system the 
directed polymers in random environment. 

Coming back to BRWRE, fix the environment q = {qt,x] {t,x) G N x Z*^} 
with mt^x > 0, Q-a.s. Set exp(/3?7t^a;) = mt,x for each {t,x) G N x Z'^. Then we 
have from Lemma 1.1 that 

E'^ [Nt,x] = E% [eMPHt - tXi(3) ):S'^ = x] and E^ [Nt] = Zt , 

where E'^[-] denotes the expectation with respect to P'' since we have that 

— ;y^%-ix'(?): if |x — vl = 1, y/x = fc, 

0, otherwise. 

Here we remark that A(/3) = log(m) so we can construct DPRE from BR- 
WRE. In [20] we find the converse, that is, how to construct i.i.d. random 
offspring distributions qt^x from the environment rjt^x- 

2. Proof of Theorem 1.3. 



2.1. Preparation. First we introduce some useful notations. We define 
u;(!«,^,y,y) for = (xi,X2),x= (xi,X2),y = (yi,y2),y = (yi,y2) G N x Vn by 
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' 0, if fly/^fly/^ = 0, 

1, if xi /xi,ay/^ay/^ 7^0, 

0, if x = i,yi /yi,ay/^ay/^ /O, 
- -1 



= < 



j>mm{k,l} 



Let {St, Pa) be an independent copy of (5t, Pg) and P^'^ be the product mea- 

sure of Pg and P| for !!^, ^ € Z"' x Vn- Then we have the following Feynmann- 
Kac formula: 

Lemma 2.1. For t e N, (y,y), (y,y) G Z'^ x Vn, 

(2.2) = ^S^"*' '^'^ = 

where et is defined by 

t 

et = JJ w{Su-i,Su-i,Su, Su)- 



u=l 

Proof. From Lemma 1.1 we can write that 

NlyNly = El[Ct; St = {y, y)]£;^[G; 5^ = {y, y)] 
= i?J°[GG;(5t,5t) = ((y,y),(y,y))], 

where C,t is Ct defined by St- It is easy to see that the P-expectation of 
the right-hand side coincides with the right-hand side of (2.2) from Fubini's 
theorem. □ 



By using this formula we can represent the uniform square integrability 
of Nt in terms of the environment, that is, {qt,x', {t, a;) E N x Z*^}. This is the 
same condition as in [20]. 

Lemma 2.2. Suppose d'>2>. Then the following are equivalent. 
(i) svLY>t>iE[{Nt)^]<^- 
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(ii) m > l^m^"^^ < ooand a = ^^^^-f^ ^ < where iTd is the return probabil- 
ity of a simple random walk in Z,'^. 



Proof, (i) =^ (ii): From Lemma 2.1 we can write that 
(2.3) Emf] = Elp,]. 

It follows from Fatou's lemma that 



liminf 



<snpE[{Nt)^]< oo. 

t>i 



By definition we can see that et = et+i on {Si / Sj} almost surely. The 
random walk — is transient since it is irreducible on Z'^ for d>3 and 
hence the limit eoo = hnit^-oo exists P^'^-almost surely. Let r be the first 
splitting time 

(2.4) T = M{t>l;St^St}. 
Then it is easy to see that 

El:'-[er;Sl^Sl] = 
since t(j(x, i, y, y) = when x = x and yi ^yi. This implies that 

(2.5) e, = 0, ^°;°-a.s. on {5^/5^}. 

If S'^ 7^ S^, then ^ S'^ ior u<t and, therefore, it is clear that 'w{St, St, St+i, 
St+i) depends only on Sj — Sj, S^_^i/Sf and S^_^i/ S"^ iov u<t {shift invari- 
ance). From this we have that 

(2.6) -^^'^600] is constant for x= (7/,X2),i= (y,X2) G Z'^ x Vn,X2 /x2- 

From (2.5), (2.6) and Markov property we can deduce that 

00 

k=l 

= T.Es>r;r = k]E^--[e^], 

k=l 

where -E"^'^ denotes the expectation with respect to P^'^ for ^ = (0,X2),^ = 
(0,X2), X2 /x2 E Vn- In the following we use E^^^- in this sense. 
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It is easy to see that E^^^~[et]T >t] = Indeed, for t = 1 we have that 
T > 1] = E^s°s^w{0, 0, Si, Si); Si = Si] 
= E^s^g[{0'Si/o)~^]Si = Si] 

= m-^ (■.■(1-6)). 

By induction we have from Markov property that 
[et; r >t] = Elf-[et; Sj = S„j = 1, . . . ,t] 

(2.8) = [et_i; 5,- = Sj,j = 1, . . . ,t - 1]^^°'° h; 5i = ^i] 
= m~K 

Also, it is easy to see that E^^^-[ei;Si / Si] = ?7T,~^(m^^) — m). Indeed we 
have that 

= E E i^°;°[«^(0,0,y,y);5i = (yi,y),5i = (yi,y)] 
= ?n~2^ E [™a^{VO'«y/o}]^"^VoOy/0 

(2.9) =m-2J^min|j^g(j),5](?(i)| 

fc^^ i>fc j>£ ' 

k>l j>k 

= m~'^{mP'^ — m). 
Prom this we can calculate E^'^~[er',T = t] as follows: 

Elpr;r = t] = Elpt-il{r>t~i}E'^'s^ei;Si + ~Si]] 

= ?Ti~^*^"^^m~^(m^^^ — m). 

Later we will prove 

(*) ^5 [eoo-- S}^ sit Gn*]=a{l-7ra)>0. 

These imply that m> 1 and m^^^ < oo from (2.7). In the remainder we check 
(*) and that E^'^Jcoo] < oo implies that a < —. 
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We divide E^g^-[eoo] according to the number of meetings of two random 
walks {Si Si). ' 

oo 

-^^'^600] =Xl-^sy^~''^^ ^ oo,r£+i =00], 

where we define tq = 0, T£ = inf{t > Ti^i^Sj = S}} for ^ > 1 with inf = 
+00. We can obtain that 

^5|[en;ri < 00] = i?°\°^,[i?5;|[e.J;ri < 00] 

= ^5i^5i[a;n <ooj 

To justify these equalities we first remark that 'w{St-i, St-i, St, St) = 1 for 
2 < t < Ti < 00 P^'^-a.s. So we can write = w{So, Sq, Si, Si) P^'^-a.s. and 

we see its P^f^^-a.s. independence of {S^,S^) from (2.1). Next we have that 

-iX2,X2r„ 1_ 77X2, X2[ 



-^52^2 [^n] - E'^l'^l[w{So, So, Si, Si)] 

, . _ -sp -^E»>fcgo,o(^)Ej>^go.o(i)] E»>fcg(^)Ej>£g(i) 

^ ■ ^ ",^1 E.>fe9(0E,>.9(j) 

(2.11) =a. 
Also we know that 

P^/^-, [n < 00] = P^;^-, [5,1 = Sl^ t>l] = PliSi = 0,^ t > 1] = vr,. 

These imply (*). Also, it follows from Markov property that 

00 

£^5;|[eoo] = ^(^jyeri ; n < cx)])^^g[eoo; Ti = (X)] 

00 

= y^(a7rd)^a(l - tt^) < 00, 

and therefore this implies that ana < 1. 

(ii) =^ (i) This has been proved in [20], Theorem 1.1, page 1623. □ 



The next theorem means the delocalization (see the remark after the 
proof). 
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Theorem 2.3. Suppose d>3 and (1.11). Then there exists a constant 
C such that 

(2.12) E^f-\et;Sl = S}] < Ct-^'"^ for all t € N. 

Remark. This theorem has been aheady proved in [20], Proposition 
1.3, page 1624, but we prove it in this article by another way because it 
contains a certain important estimate which is used in the proof of our main 
theorem. 

Proof of Theorem 2.3. From the same argument as in the proof of 
Lemma 2.2, we can obtain that 

t 

£;°;°[ei; si = S}\ = Elpt; r>t] + Y, r = k, Sj = S}] 

k=l 



m * 



k=l 



where c is the constant given by m ^(m^^^ — m) and where in the last term 
we used (2.8) and (2.9). It is clear that m~* + m~*+^c < Ct~^^'^ and hence 
it is enough to estimate E^^^~[et-k] Sj_j^ = Sj_f,]. By using tj, j > 0, we can 
rewrite it as 

t-k 



s,s 
e=i 

t-k 

= E''^''g[et~k;Tl=tl,T2-Ti=t2,...,Ti-Ti_i=ti]. 
i=\ yt^=t-k 

If we set at = i?^'^ [e^-i ; ri = t], then it follows from Markov property and 
shift invariance that 

E^^^et-.k\T\ =ti,T2-Ti = t2,...,T£- rf„i = ti] = at^at^ ■■■ati, 

if ti H \-ti = t — k. We remark that from (2.10) 

at = aPlf^^{n = t)<c,t~''/^ 

(2.13) '^0't = -^^y^ri ; Ti < oo] = 7? = aiTd < 1 and 



t>i 



Y at-,---at, = (£^J^ [cr^ ; Ti < oo] )^ = 7/^ 
t>iti+---+ti=t 
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where we used on the first fine the fact that sup2.g2d Pgi [Sj = x] = 0{t~'^^^). 
From these properties we prove that there exist /3 < 1 and Ci > such that 

(2.14) at,--- at, < Cil3H~'^/'^ for aU t > 1. 

We consider the sequence {cfc}fc>i satisfying that for < e < 1, 
(2-15) Cfc+i = (13^^' + 

where ci is given in (2.13). First we will prove for all > 1 the following 
inequality holds: 

(2.16) at,--- at^ < Ckt'"^^^ for ah t > 1. 

ti+-+tk=t 

Indeed this inequality holds for k = 1. Suppose (2.16) holds for A; > 1. Then 
we have the following inequality from (2.13): 

«ti • • • atk+i 

ti+---+ifc+i=t 
t-i 



J2( Yl ah---atAat- 

o— V, J l~t, =8 / 



s=k HiH Hk=s 

Y1 at,---at\i{t-s)-'^/^+ Y CkS-'^/^at-s 

s<et Hk=s ^ et<s<t 

E at,---atX,{t-et)-''/'+ Ckiety'^t-s 

S<£t\l-\ Hk=S ^ £t<S<t 

and hence (2.16) holds for k + 1. We choose e such that rj < e'^l'^ < 1. Then 
we have < C(^^)^ for all /c > 1 by simple calculation and (2.14) follows. 
Therefore, we obtain that 



Eg[e,_,; 5i_, = Si,] < J] C,/3'{t - ky^l^ <C2{t- ky 

and from this it is easy to check (2.12). □ 

Remark. We define pt and TZt by 

pI = msixpt(x) and TZt = Pti^)- 



■d/2 
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is the density at the most populated site while TZt is the probability that a 
given pair of particles at time t are at the same site. Clearly (p^ )^ <T^t ^ Pt- 
The above theorem can be interpreted as if we suppose that d>3 and (1.11), 
then 

TZt = 0{r'^l'^) in P(-|iVoo > 0)-probability. 
This can be seen as follows: 

and linit-j-oo-^t = -^oo > 0, P{-\NaQ > 0)-a.s. However, we know from Lem- 
ma 2.1 and Theorem 2.3 that 



E E 

and hence we have that 



E'f^[euSl=~Sl] = 0{t-'^'^) 



E 



since P(iVoo > 0) > 0. 



2.2. Some propositions. We now show Theorem 1.3 by using the argu- 
ment in [4]. First we introduce some notations. Let {^t}t>i be i.i.d. random 
variables with values in W^. We denote by Xt a random walk whose steps 
are given by the ^^'s. Moreover, we assume that -E[exp(0 • ,^i)] < oo for 9 in 
a neighborhood of in R'^. We define p{9) by 

(2.17) p(^) = lni?[exp(0-6)]. 
Then it is obvious that 

exp(0 • Xt - tp{e)) 

is a martingale with respect to the filtration of the random walk. 

We win use standard notation x" = • • • x""/ and (^)" = {-^T^ ■ ■ ■ (^) 

for n = (ni, . . . , n^) G N'^ and x G M*^. For n = (ni, . . . , Ud) G N'^ the polyno- 
mial Wn{t,x) is defined by 

d_ 

where |n| = ni + ■ ■ ■ + rid- We write 

(2.18) Wr,{t,x)= Yl ^n(i,i)xV. 

(i,j)GN'*-ixN 

The coefficients An{i,j) depend on the derivatives of p in 0. The following 
lemma gives some useful properties of Wn{t,x). 



Wnit,x)=[ — ] expie-x-tp{e))\,^o, 
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Lemma 2.4. For a general random walk with exp{p{6)) < oo for 6 in a 
neighborhood of and E[S^i] = 0, we have: 

(a) +2i > |n|, then An{i,j) = 0. 

(b) Coefficients with |i| + 2j = |n| depend only on the second derivatives of 
p at 0, that is, on the covariance of ^i. 

(c) //|i| = |n|, i/ien An(i,0) = (5ii,„i(5j2,„2 • • -Jj^^n^. 

Proof. We have that {^){x ■ 9 - tp{6))\g=o = Xi and {-§gY{x ■ 6 - 
tp{e))\e=o = -*(M)'(p(^))k=o for |i| > 2 since ^^p{9)\9=o = 0. (a)-(c) follow 
from Faa di Bruno's formula [8], Theorem 2.1, page 505, and from the fact 
that ^e^U=o = 1 for all A: G N. □ 

Wn{t,Xt) is a martingale with respect to the filtration of the random 
walk. Coming back to Markov chain (S,Pg) we have that 

(2.19) Y^it) = E''s[Wn{t,Sl)Ct] 

is an J-j-martingale since Ct is an "Ht-martingale. Indeed we have that for 
any set B £ Ft-i, 

EA[El[W^{t, SDCt] ■■ B] = E'siWnit, Sl)E[Ct : B]] 

= E^s[Wr.it,Sl)E[Ct^,:B]] 

= EA[E's[W^{t,Sl)Ct^i]:B] 

= EA[E's[W^{t-l,Sl,)Ct^i]:B]. 



Proposition 2.5. Suppose d>3 and (1.11). Then we have that for 
each n G N"^ with |n| / 

lim t~l"l/2y„(t) = 0, p-a.s. 

t— >oo 

Proof. We show that the J^^-martingale 

Z,'^j2s~\-\'\Y^{s)-Y^{s-l)) 

s=l 

remains L^-bounded. This implies that Zt converges a.s. and hence Proposi- 
tion 2.5 follows from Kronecker's lemma for |n| ^ 0. For simplicity we write 
Wr,{t,S}) = Wr,{t,S). It is enough to show that E[{Yr,{t) - Y^{t - < 
(j^\n\-d/2_ injiee^i obvious that 

oo 

(2.20) snpE[Z^] = V s-l"li^[(y„(s) - Yn{s - 
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and hence, if we show that E[{Yn{t) - Ya{t - < Ctl"!"'^/^, then the 
right-hand side of (2.20) is finite. We can write that 



(2.21) 



E[(,Yn{t) - Yn{t - l)f] = E[{E''s[Wn{t, S)Ct - Wr,{t - l,S)Ct~i]f] 

= E[iE''s[WUt,S){Ct-Ct-i)] 

+ E^siiWnit, S) - W^it - 1, S))Ct-i]y 
= E[{E^s[Wn{t,S){Ct-Ct^i)])\ 



where we use the fact that E^[{Wn{t, S) -Wn{t- 1, 5))Ct-i] = P-a.s. from 
the observation after the proof of Lemma 2.4. Moreover, we have that 

the right-hand side of (2.21) 



£;[Eg[W^„(t,5)(Ci-G-i)]i?a[W^n(i,5)(Ct-G_i)]] 



E 



0,0 



(2.22) 



W^it,S)Wnit,S) 



X E 



Ct-iCt 



\^St/St 



A 



t,St. 



''St/St 



= E''^-^[W^{t,S)W^{trS){et-et.MSU=~SU}] 

= El''-[\W^{t, S)\\et - et-.i)l{Sl, = Sl,}], 

where we used on the last line the following facts obtained from Markov 
property and (2.1): 



E 



A St 



/I St 



^St/St-1 



et-i{w{St-i, St-i, St, St) - 1) 



1 ^^'g-a-s. on {S}^^ / S^^^}. 



Kast/St-i 
and w{St^i,St^i,St,St) 

It is easy to see that \Wn{t,S)\'^ < CglS^Lipl"! + Qtl"! from Lemma 2.4, 
where C3 and C4 are constants dependent only on n and d. We have already 
proved that £;°'^[et;5t^ = S}] < Ct''^''^. Therefore, from (2.22) we have to 

estimate the values E^f-^[\Sl\^\^\etl{Sl = S}}] and E^f^[\Sl\'^\'^\ et+il{Sl = 

S}^]. However, we know from Markov property that 

S]}] 



E'f~,[\Sl\'\-\et^,l{Sl 



S,S^ 



< max 



m 



m 



s,s 

2 ,a\El'l[\Sl\'^'^\etl{Sl 



(2) 



w{St, St, St+i, St+i)]] 
Si}], 
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where we used the fact that 



[a, ify = y,y/y. 

Therefore, it is enough to show that E^fj\Sl\^\'^\ etl{Sl = S}}] < Ct\^\~'^/'^. 
We define Uk for A; G N by 

o"o = inf{t > 0; St 7^ St} and (7^ = inf{t > cr^.i; S/ = Sj } for /c > 1 

with inf = +00. We remark that (Tq = t where r is defined by (2.4). More- 
over, let Xto,ti,-,tk = l{^o = tQ, fJi - do = ti, CTfc - ak-i = tk}. Then with 
a similar argument to the proof of Theorem 2.3, we can write that 

[l^tT'-'et; Si = S}] = E'^-[\Sl\'\-W,oo > t] 

k=Oto + ---+tk=t 

Since jSt^pl"! < it is clear that 

(2.23) i?°'°[|S,Yl"let;ao>t]<t2'"l/"i* 

and hence we have that [jSt^^l^let; ao > t] < Ctl"!-'^/^ 

In the remainder we will show that there exists a certain constant C > 
such that 

fc=Oio+---+tfe=t 

If this has been shown then we complete the proof of Proposition 2.5. Since 

k 
1=1 

it is obvious that 

(2.24) Y: E ^55[l^*l""'^*X*o,...,tJ 
k=oto+---+tk=t 

t k 

(2.25) <5:(^+i)""'E E ^5:°[i^^,-5LJ""'«*x*o,...,.j 

k=0 i=l to + ---+tk=t 

t 

(2.26) +^(fc+l)2|n| ^ i?°;°[|5^ol""'e*Xto,...,*J- 
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By using Markov property and shift invariance we have that, for 1 <i <k, 



foH \-tk=t 



u/, —t ' ^^-^n ' / 



(2.27) 

xii;^'J[|5ij2He.,;ai = t,] 



...4-/, =f ^^'^n 4 / ' 



< 

to+---+ifc=t "-j^OA 

where c is the constant given by c = m~'^(m^'^^ — in). It is easily seen from 
(2.10) that 

^3[|5/pl"lei;ai =t] =i?j|[|5i|2|"l;ai =t]a 

where we have used on the third hne the fact that sup^.g^d Pgi [S] = x\ = 
0{t-'^/'^) [19]. Therefore, it follows that 

the right-hand side of (2.27) 

to+U<t 

where we use (2.14) and the fact that m~* < Ct~'^/'^ . Since this inequality is 
independent of z, we have that 

oo 

the right-hand side of (2.25) < E<^('^ + l)2|"l+i/3^-iil"l-^/2 

k=l 
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where C is a constant depending only on n and d. A similar argument holds 
for the right-hand side of (2.26). Indeed we have that 

io+-+tfe=t 

to+-+tk=t ' j^O 

< E cm-(*o-i)tflj]a,^, 

to+-+tk=t jVO 

to<t 

where we use (2.14) and the fact that t^l^'/jn* < Cgt"'^/^. Hence we can 
obtain that 

the right-hand side of (2.26) < Clo^"''/^ 

where Cio is a constant depending only on n and d. From these we have 
that 

the left-hand side of (2.24) < Ctl"!-'^/^^ 

so that 

k=Oto+---+tk=t 

where C is a constant depending only on n and d. Hence the proof is com- 
plete. □ 

Since we have proved Proposition 2.5 we can show Theorem 1.3. 

2.3. Proof of the result. From [14], Theorem 3, page 363, it is enough to 
show the following proposition instead of Theorem 1.3. 

Proposition 2.6. Suppose d> 3 and (1.11). Then for all n = {ni, ... ,nd) £ 



(2.28) lim E^s 



1 d =Noo [ x"di/(x), p-a.s., 

J J JRd 



where v is the Gaussian measure with mean and covariance matrix 
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Proof. By induction it follows from Lemma 2.4(a), (c) and Proposi- 
tion 2.5 that for all n G N'^ 



(2.29) 



sup 

t>i 



Eg 



Ct 



< oo, P-a.s. 



To see this we divide Yn{t) into three parts as follows: 

y„i(t) = ti"i/2E0 



(2.30) y„2(^)^^|n|/2^0 



|i|+2i=|n|,j>l - 



5, 



1 \ 1 



and 



Then we can write 



1 \ ' 1 

Ct 



(2.31) 



|i|+2i<|n 



Ct 



^-\n\/2Yl^t) 



= t-\"\/^{Yn-Y^-Y^). 

We suppose that (2.29) holds for n G N'^ with |n| < k. From Proposition 2.5 
we have sup^>i Yn(i)| < oo P-a.s. for all n G N"^. It is easy to check 

that for n G N*^ with |n| = A; + 1, 

supt-l"l/2|y„2(t)| <oo and supt~\''\/^\Y^(t)\ < oo, P-a.s. 
t>i t>i 

Thus (2.29) holds for ah nG N*^. Therefore, we conclude that 

(2.32) lim rl"l/2y„3(t) = 0, P-a.s., 

and hence from (2.31) and Proposition 2.5 that for |n| > 1 

(2.33) lim t~l"l/2(y^i(t) + Y^{t)) = 0, P-a.s. 

On the other hand, let Z be an M'^-valued random variable with density u. 
Then it can be seen that pi(^) is a polynomial of degree 2 where pi{0) is 
given by (2.17) for = Z. Moreover, we have that for |n| > 1, 



E 



|i|+2i<|nl 



ALMOST SURE CLT FOR BRANCHING RW IN RANDOM ENVIRONMENT 23 



where A'^{i,j) is defined by (2.18). From Lemma 2.4, A'^{i,j) corresponds 
with An(i,i) for (i,j) with |i| + 2j = |n| and hence we can write for |n| > 1 



(2.34) 



E 



|i|+2i=|n|J>l 



0. 



Here we remark that A'^{i,j) = for (i,j) with |i| + 2j < |n| since {-^Y x 
pi{e)\e=o = for j G N'^ with |j| > 3._ 

We know that limt^oo Eg[(t] = ^oo for |n| = which gives (2.28) for 
|n| = 0. If (2.28) holds for ah n G N'' with |n| < k, then we have that for ah 
n e N*^ with |n| = A; + 1, 



(2.35) hm t-l"l/2y^2(t) = NooE 



|i|+2i = |n|, 
J>1 



P-a.s. 



From this, (2.32) and Proposition 2.5 it fohows that the right-hand side of 
(2.31) converges to 



-Nr^E 



|i|+2i = |n|, 

i>i 



almost surely as t oo, so that (2.28) holds for n G N'^ with |n| = k+1 from 
(2.34). Therefore, we complete the proof of Proposition 2.6 and Theorem 1.3. 
□ 
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